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We put forward a scenario where blazars are classified as flat-spectrum radio quasars, BL Lacs, low synchrotron, 
or high synchrotron peaked objects according to a varying combination of Doppler boosted radiation from the 
jet, emission from the accretion disk, the broad line region, and light from the host galaxy. We thoroughly test 
this new approach, which builds upon unified schemes, using Monte Carlo simulations and show that it can 
provide simple answers to a number of long-standing open issues. We also demonstrate that selection effects 
play a very important role in the diversity observed in radio and X-ray samples and in the correlation between 
luminosity and the peak frequency of the synchrotron power (the so-called "blazar sequence"). It turns out 
that sources so far classified as BL Lacs on the basis of their observed weak, or undetectable, emission lines 
are of two physically different classes: intrinsically weak-lined objects, more common in X-ray selected samples, 
and heavily diluted broad-lined sources, more frequent in radio selected samples, which explains some of the 
confusion in the literature. 



1. INTRODUCTION 



Blazars are radio loud active galactic nuclei (AGN) 
pointing their jets in the direction of the observer 
[1, |48[ . While all blazars emit variable, non-thermal 
radiation across the entire electromagnetic spectrum, 
they come in two main subclasses, whose major differ- 
ence is in their optical properties: 1) Flat Spectrum 
Radio Quasars (FSRQs), which show strong, broad 
emission lines in their optical spectrum, just like radio 
quiet quasars; and 2) BL Lacertae objects (BL Lacs), 
which are instead characterized by an optical spec- 
trum, which at most shows weak emission lines, some- 
times displays absorption features, and in some cases 
can be completely featureless. Historically, the sepa- 
ration between BL Lacs and FSRQs has been made 
at the (rather arbitrary) rest-frame equivalent width 
(EW) of 5 A [3, IHJ]. However, no evidence for a bi- 
modal distribution in the EW of the broad lines of 
radio quasars has ever been found and, on the con- 
trary, radio-selected BL Lacs were shown to be, from 
the point of view of the emission line properties, very 
similar to FSRQs but with a stronger continuum [40| • 
Most BL Lacs selected in the X-ray band, on the other 
hand, had very weak, if any, emission lines, and, when 
studying the properties of the X-ray selected Einstein 
Medium Sensitivity Survey (EMSS) sample, another 
criterion to identify them, this time to sep arate them 
from galaxies, had to be introduced [45|]. This was 
based on the Ca H&K break, a stellar absorption fea- 
ture typically found in the spectra of elliptical galax- 
ies. Given that its value in non-active ellipticals is 
~ 50%, a maximum value of 25% was chosen to ensure 
the presence of a substantial non-thermal continuum 
superposed to the host galaxy spectrum (45{. This 



was later revised to 40% 



Blazar classification depends then on the details of 
their appearance in the optical band where they emit 
a mix of three types of radiation: 1) a non-thermal, 
jet-related, component; 2) thermal radiation coming 
from the accretion onto the supermassive black hole 
and from the broad line region (at least in most radio- 
selected sources); 3) light from the host (giant ellip- 
tical) galaxy. Figure [1] represents these three com- 
ponents as red, blue and orange lines, overlaid to 
the spectral energy distribution (SED) of four well- 
known blazars [2fJ. The strong non-thermal radia- 
tion, the only one that spans the entire electromag- 
netic spectrum, is composed of two basic parts form- 
ing two broad humps, the low-energy one attributed to 
synchrotron radiation, and the high-energy one, usu- 
ally thought to be due to inverse Compton radiation 
[e.g. The peak of the synchrotron hump (i/ pea k) 
can occur at different frequencies, ranging from about 
- 3 x 10 12 Hz to over 10 18 Hz (see 3C 273 or 3C 279 
and MKN 501 in Fig. [1} reflecting the maximum en- 
ergy at which particles can be accelerated [e.g. l2(il |. 
Blazars with rest-frame ^ pca k < 10 14 Hz are called 
Low Synchrotron Peaked (LSP) sources, while those 
with 10 14 Hz < v pcak < 10 15 Hz and u peak > 10 15 Hz 
are called Intermediate and High Synchrotron Peaked 
(ISP and HSP) sources, respectively @. These def- 
initions extend the original division of BL Lacs into 
LBL and HBL sources [30(. An expanded version of 
the work presented here is given by [181 ] . 
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Figure 1: The SEDs of four representative blazars: two FSRQs, 3C 273 and 3C 279, and two BL Lacs, MKN 501 and 
BL Lac. The lines in color denote the three main components of blazars SEDs, namely non-thermal radiation from the 
.jet (red), emission from the disk and from the broad line region represented by the composite QSO optical spectrum of 
[49l ] (blue), and light from the host galaxy, represented by the giant elliptical template of [24| (orange). The two 
vertical lines indicate the optical observing window (3800 — 8000 A). 



2. CURRENT STATUS 

The two main blazar subclasses have many differ- 
ences, which include: 

1. different optical spectra (by definition). There 
are, however, a number of BL Lac - FSRQ tran- 
sition objects, which include even BL Lacertae 
itself, the prototype of the class, which displays 
at times moderately strong, broad lines [e.g. |50| 
and 3C 279, a well-studied FSRQ, which can 
appear nearly featureless in a bright state j34|; 

2. different extended radio powers. Most BL Lacs 
have extended radio powers and morphologies 
consistent with those of Fanaroff-Riley (FR) 
type L while basically all FSRQs are FR II- 
like [48|, and references therein]. However, some 
radio-selected BL Lacs are known to posses an 
FR II- like structure [e.g. HI; 

3. very different redshift distributions. FSRQs, 
similarly to radio quiet QSOs, are typically 
found at redshifts ~ 1 — 2, and up to ~ 5.5, 
while BL Lacs are usually much closer with 
very few cases at z ^ 0.6 [23]. A large fraction 



of BL Lacs, however, [~ 43% in BZCAT and 
> 50 - 60% of the BL Lacs in the Fermi 1 and 2 
year AGN catalogs: Q], 0] have no measured red- 
shift, due to the lack of any detectable feature in 
their optical spectrum, despite the use of 8/10- 
m class optical telescopes for the spectroscopy 
identification campaign; 

4. different cosmological evolutions. Radio and X- 
ray selected samples have shown that the two 
blazar subclasses have very different cosmolog- 
ical evolutions, with FSRQs evolving strongly 
(again similarly to radio quiet QSOs) and BL 
Lacs evolving at a similar, or perhaps lower rate, 
in the radio band, or even showing no or nega- 
tive evolution in the X-ray band [e.g. Hi], H3) Ej] ; 

5. widely different mix of FSRQs and BL Lacs in 
radio and X-ray selected samples, with the latter 
typically including a much larger fraction of BL 
Lacs than the former. In fact, while only ~ 15% 
of WMAP5 blazars are BL Lacs (Section l3~Tj) . 
this fraction is instead ~ 70% in the EMSS com- 
plete sample, which includes 41 BL Lacs and 15 
FSRQs ||[33; 
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6. widely different distributions of the synchrotron 
peak energy ^ pea k- The rest-frame Vp Ca k distribu- 
tion of FSRQs is strongly peaked at low energies 
(( ^pcak) = 10 131±01 Hz) and never reaches very 
high values (^ pca k ~ 10 14 ' 5 Hz) independently of 
the selection method [13], while the t'pcak dis- 
tribution of BL Lacs is shifted to higher values 
by at least one order of magnitude. It can also 
reach values as high as f poa k ^ 10 18 Hz and its 
shape varies strongly depending on the selection 
band [that is radio, X-ray or T-rav: I2I [20j . 

Some of these differences have been explained by 
so-called unified schemes, which posit that BL Lacs 
and FSRQs are simply FR I and FR II radio galaxies 
with their jets forming a small angle with respect to 
the line of sight Radio galaxies would then be 

the "parent" population of blazars. Due to relativis- 
tic beaming this has enormous effects on their appar- 
ent emitted power and luminosity functions (LFs) and 
can explain their different extended radio powers and, 
partly, their cosmological evolutions. However, unified 
schemes per se cannot account for transition objects, 
the different evolution of radio and X-ray selected BL 
Lacs, and ^ pea k distributions. 

3. A NEW SCENARIO FOR BLAZARS 

In our new scenario, the observed blazar optical 
spectrum is the result of a combination of an in- 
trinsic EW distribution and the effects of three com- 
ponents: a non-thermal, jet related one, a thermal 
one due to the accretion disk, and emission from the 
host galaxy. Different mixes of these components de- 
termine the appearance of the optical spectrum and 
therefore the classification of sources in FSRQs (dom- 
inated by strong lines), BL Lacs (with diluted, weak 
lines, if a standard accretion disk is present), and 
radio-galaxies (where the host galaxy swamps both 
the thermal and non-thermal nuclear emission present 
in blazars) . The other novel component is a single LF 
whose evolution depends on radio power. 

3.1. Simulation ingredients 

Our idea was tested through extensive Monte Carlo 
simulations, which include the following ingredients, 
which we kept as simple as possible and tied as much 
as possible to observational data: 

1. Luminosity function We derive the LF and 
evolution of blazars at 41 GHz from the Wilkin- 
son Microwave Anisotropy Probe (WMAP5) 
sample [52| . which is practically equivalent to 
a radio-selected sample. We extend on the 
work of [T(| and define a flux-limited sample of 
high Galactic latitude sources {f^iGHz > 0.9 Jy, 



1 611 1 > 15°) including 161 FSRQs, 29 BL Lacs, 
and 10 blazars of unknown type. By applying a 
maximum likelihood technique to the WMAP5 
blazars we obtain, together with the evolution 
discussed below, a best-fit local LF $(P) oc P~ 3 
(in units of Gpc~ 3 P -1 ) between 1.9 x 10 24 and 
4.2 x 10 27 W/Hz, which we assume in our radio 
simulations. Throughout this paper we use a 
ACDM cosmology with P = 70 km s _1 Mpc~\ 
tt m = 0.27 and fl A = 0.73. 

2. Cosmological evolution Powerful (P r ^ 10 26 
W/Hz) radio sources display strong evolution at 
low redshifts followed by a decline at higher red- 
shifts [e.g. [5l|. We parametrize this behavior 
with a model of the type P(z) = (1 + z) k+l3z , 
which allows for a maximum in the luminosity 
evolution followed by a decline. A maximum 
likelihood technique applied to the WMAP5 
sample allows us to derive k = 7.3 and (3 = —1.5 
in the — 3.4 redshift range (which implies a 
peak at z ~ 1.85), which we assume in our sim- 
ulations. Lower luminosity (P r ^ 10 26 W/Hz), 
mostly FR I radio sources display a much weaker 
cosmological evolution, which reaches « zero 
at P r < 10 25 W/Hz [e.g. [H, and references 
therein]. We took this into account by using 
the radio LFs of BL Lacs and FSRQs derived 
from those of FR Is and FR lis and based on 
the beaming model of (48|. We then used the 
fraction of beamed FR I blazars in bins of radio 
power to simulate the fraction of non-evolving 
radio sources as a function of power. 

3. Non-thermal component To represent the 
non-thermal/jet component, we assume a simple 
homogeneous synchrotron self-Compton model 
[SSC, see, e.g. |47|, and references therein] with 
relativistic electrons distributed as a power law 
at low energies and as a log-parabola at high en- 
ergies [1^. This model represents well the syn- 
chrotron part of the observed SEDs, which al- 
ways extends at least to the optical band where 
the classification of blazars as FSRQs or BL 
Lacs occurs. As for the inverse Compton emis- 
sion, which can be important in the soft X-ray 
band, we set the Compton dominance so as to 
reproduce the observed f x /f z in FSRQs. The 
Lorentz factors of the electrons radiating at the 
peak of the synchrotron SED component (7 pea k) 
cover the range ~ 10 25 — 10 4 5 , which is that ex- 
pected for typical parameters of the SSC model 
as shown in Fig. 36 of [2j. The shape of the dis- 
tribution was chosen to reproduce the observed 
^pcak distributions in radio and X-ray selected 
samples of blazars. We assumed a mean value of 
15 for the Doppler factor S. This was chosen to 
be consistent with the mean superluminal speed 
Papp ~ 12 obtained by (HI, and with the typical 
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Lorentz factor V ~ 15 derived by [2l[ (since for 
the angle that maximizes the apparent velocity 

s ~ /3a PP ~ r). 

4. Accretion Disk and broad emission lines 

We use the quasar spectral template of [4j| [see 
Fig. 4 of [20|. A standard accretion disk is 
likely to be present only in so-called "high ex- 
citation" radio galaxies (HERGs), while it ap- 
pears not to be there, or be less efficient, in 
low-excitation ones (LERGs). Almost all FR 
Is are LERGs, while most FR lis are HERGs, 
although there is a population of FR II LERGs 
as well. Observational evidence suggests that 
LERGs, which means all FR Is and some FR 
lis, either do not possess an accretion disk, or if 
the disk is present is much less efficient than in 
FR lis (i.e. of the Advection Dominated Accre- 
tion Flow [ADAF] type) [e.g. SI . We have then 
associated the presence of a standard accretion 
disk only with beamed FR II sources and as- 
sumed that all those with an FR I parent (the 
non-evolving sources) have no disk. 

5. Equivalent width distributions The in- 
trinsic (before dilution) distributions of the 
EW of the broad lines (Lya, C IV, C III, 
Mg II, H/3, Ha) have been assumed to be 
those of a sub-sample of the radio quiet 
QSOs included in the SP SS DR7 database 
(| http:/ /www. sdss.org/dr77| . We assumed 
Gaussian distributions characterized by the 
measured means ((EW Ho )= 200 A, (EW H ^)= 
23 A, <EW Mg -ii>= 18 A, (EW C -iii)= 16 A, 
(EW C -iv}= 20 A, (EW Ly a)= 47 A) and dis- 
persions for the various lines. 

6. The disk to jet power ratio The disk and jet 
components in blazars are known to be corre- 
lated and possibly of the same order of magni- 
tude [see, e.g. ahj. We are interested in the 
somewhat simpler question of determining how 
the luminosity of the accretion disk (blue bump 
intensity at 5000 A) scales with radio power (at 
5 GHz). The relevant data were derived by us- 
ing the very large amount of multi-frequency in- 
formation included in public databases and the 
tools that are now available to analyze SEDs 
,46j . Figure [T] gives some examples of represen- 
tative objects. The amount of thermal flux in 
each FSRQ was estimated by matching the com- 
posite optical QSO spectrum of to the SED 
data in the blazar rest-frame. Fig. [1] gives ex- 
amples of the matching of the composite QSO 
spectrum (blue line) to the data for the case of 
3C 273 or 3C 279. An upper limit was esti- 
mated for BL Lac objects by placing the com- 
posite QSO spectrum in the SED at an inten- 
sity such that the optical lines would not be 



detectable in the optical spectrum (typically a 
factor ten below the observed flux. Fig. [T] il- 
lustrates the case of MKN501). This was done 
for the blazars selected in two surveys, WMAP5 
and the EMSS. We found that otr—BiueBump (the 
slope between the 5 GHz luminosity and the blue 
bump intensity at 5000 A and defined by Ldisk = 
LAv ; /^5Gffz)~ Qr ~ BIUcBump ) correlates with 

\ 5000A'_ °^ nz > I _ 

radio luminosity, although with a large scat- 
ter. We then fitted a simple linear relation- 
ship between the two variables {ct r -BiueBump = 
0.04*log(L rac iio) — 0.39) and assumed a Gaussian 
distribution around it with a dispersion of 0.1. 

7. Host galaxy We have assumed that the host 
galaxy of blazars is a giant elliptical with fixed 
absolute magnitude of M R = -22.9 dHHO]. For 
the spectral shape we used the galaxy template 
of [24j |. who derived it combining the data from 
28 local elliptical galaxies observed in the wave- 
length range 0.12 — 2.4 /im. Figure Q] shows this 
template superposed to the SED of four well- 
known blazars. 

3.2. Simulation steps 

Our Monte Carlo simulations start by drawing a 
random value for the radio luminosity and redshift 
based on the luminosity function and evolution de- 
scribed above and a value of the Lorentz factor of 
the electron radiating at the peak of the synchrotron 
power (7 P eak) from the assumed distribution. We 
then calculated the peak of the synchrotron power 
in the source rest frame by assuming a simple SSC 
model (vpeak = 3.2 x 10 6 7p Cak £W). The magnetic field 
was fixed to B = 0.15 Gauss and the Doppler factor 
5 was randomly drawn from a gaussian distribution 
with (S) = 15 and er= 2. We then calculated the 
observed radio flux density (from the radio luminos- 
ity and redshift) and the non-thermal emission in the 
optical and X-ray bands under the assumption that 
the spectral shape of the observed emission is a log 
parabola around i/peak and that the low energy part 
of the SED (cm and mm wavelengths) is a power law 
(as is typically seen in blazars: see e.g. Fig. [T] for 
some representative examples). We then added an 
accretion (blue bump) component as described above 
(only for beamed FR II sources), re-scaling the SDSS 
quasar template to this value and drew a random value 
of the equivalent width of Lya, C IV, C III, Mg II, 
H/3, Ha starting from the EW distribution observed 
in the SDSS radio quiet QSOs. The host galaxy op- 
tical light was added assuming a standard giant el- 
liptical and we then calculated the total optical light 
and the observed equivalent width of all the broad 
lines considered by taking into account the dilution 
due to the non-thermal and host galaxy optical light. 
Sources were classified as FSRQs if the rest-frame EW 
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of at least one of the broad lines that enter the opti- 
cal band in the observer frame (which we assume to 
cover the 3, 800 - 8, 000 A range) was > 5 A. Other- 
wise, the object was classified as a BL Lac, unless the 
host galaxy dominated the optical light causing the 
Ca H&K break to be larger than 0.4 [I2,[25|], in which 
case the source was classified as a radio galaxy. A BL 
Lac whose maximum EW is < 2 A, or for which the 
non-thermal light was at least a factor 10 larger then 
that of the host galaxy 35[, was deemed to have a 



Table I Results from a simulation of a radio flux density 
limited survey (0.9 Jy) 



redshift which cannot be typically measured. 

It is important to stress that the scope of our sim- 
ulations is not to reproduce all the observational de- 
tails. While that could be possible in theory, in prac- 
tice it would require a large number of parameters and 
some speculations. Our approach is instead to keep 
the number of assumptions to a minimum, to obtain 
robust, almost model-independent conclusions. 

We are confident that our main results are stable 
to changes in input parameters. The adopted LF and 
evolution were in fact varied by la from the WMAP5 
best fit and we also used as LF the sum of the BL 
Lac and FSRQ LFs based on the beaming model of 
. We also run the simulations with values of (S) in 
the 5 to 20 range, and considering also a dependence 
of 5 on radio power [e.g. |21| . No major changes were 
obtained as compared to our default assumptions. 

3.3. Simulations of radio and X-ray 
surveys 

We simulated a radio flux density limited survey 
with / > 0.9 Jy, to match the WMAP5 sample, 
and an X-ray flux limited survey down to 5 x 10~ 13 
erg cm -2 s _1 in the 0.3 — 3.5 keV band, in order to 
be able to compare it with the EMSS. To ensure good 
statistics each simulation run included 10,000 sources. 
In the X-ray case, since radio powers reach lower val- 
ues than in the radio case, we extrapolated the radio 
LF down to 1.9 x 10 23 W/Hz assuming the same slope. 



4. COMPARING SIMULATIONS AND 
REAL DATA 

4.1. Radio flux density limited survey 

Table U summarizes our main results by giving the 
number of sources per class, their mean redshift, and 
(V/V m ) (where V is the volume out to the source 
and Kn is the volume at the distance where the ob- 
ject would be at the flux limit: [HI)- The number 
in parenthesis refers to the BL Lacs with measurable 
redshift, to which the mean redshift and (V/Vm) per- 
tain. About 3/4 of our sources are classified as FSRQs, 
with the fraction of BL Lacs being ~ 19.8% of blazars, 
which is consistent with the value of 15.3^3 q% in the 



Source type 


Number 




(V/V m ) 


FSRQs 


7,587 


1.24 


0.64 


BL Lacs 


1,879 (1,191) 0.87 


0.60 


Radio galaxies 


534 


0.04 


0.48 


Total 


10,000 


1.13 


0.63 



o 

CO 




Figure 2: Top panel: the redshift distribution of the 
WMAP5 FSRQs (solid histogram) compared to that of 
FSRQs in a simulation of a radio survey (dashed 
histogram). Bottom panel: the redshift distribution of 
the WMAP5 BL Lacs (solid histogram) compared to that 
of BL Lacs in a simulation of a radio survey (dashed 
histogram) 



WMAP5 sample. A small fraction (5%) of the sim- 
ulated blazars are classified as radio galaxies. These 
are bona-fide blazars misclassified because their non- 
thermal radiation is not strong enough to dilute the 
host galaxy component. The mean redshift for our 
simulated FSRQs agrees with the WMAP5 value of 
1.13, while for BL Lacs this is slightly larger than 
the WMAP5 mean (0.55). Fig. [2] shows the overall 
good agreement between our simulated redshift dis- 
tributions (where we have only included sources with 
a measurable redshift) and the observed ones. 63% 
of our BL Lacs have a redshift determination, in ex- 
cellent agreement with the WMAP5 value of 69^20%. 
79% of the BL Lacs (68% of those with redshift) have 
a standard accretion disk and are therefore broad- 
lined but are classified as BL Lacs only because their 
observable emission lines are swamped by the non- 
thermal continuum. As regards (V/V m ), our simu- 
lated mean values agree with the observed ones of 
0.62 ± 0.02 and 0.63 ± 0.05 for WMAP5 FSRQs and 
BL Lacs respectively [cf. also the value of 0.60 ± 0.05 
for the 1 Jy BL Lac sample: Hij]. Fig. |3] compares 
the distributions of t'peak, the synchrotron peak cn- 
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Figure 3: Top panel: the fpeak distribution of radio 
selected FSRQs taken from the work of 20] (solid 
histogram) compared to that of FSRQs in a simulation of 
a radio survey (dashed histogram). Bottom panel: the 
i/peak distribution of the BL Lacs in the radio sample of 
[201 ] (solid histogram) compared to that of BL Lacs in a 
simulation of a radio survey (dashed histogram). 
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Figure 4: Top panel: the i/peak distribution of the X-ray 
selected FSRQs in [2p| ] (solid histogram) compared to 
that of FSRQs in a simulation of an X-ray survey 
(dashed histogram). Bottom panel: the J/ pea k distribution 
of the BL Lacs in the X-ray flux limited sample of [2(| 
(solid histogram) compared to that of BL Lacs in a 
simulation of an X-ray survey (dashed histogram). 



Table II Results from a simulation of an X-ray flux 
limited survey (5 x 10 -13 erg cm -2 s" 1 ) 

Source type Number (z) (V/Vm) 

FSRQs 2,836 1.23 0.65 

BL Lacs (all) 5,622 (4,460) 0.36 0.51 

BL Lacs (log ^ pcak > 16.5) 927 (895) 0.33 0.45 
BL Lacs (log v peak > 17) 185 (177) 0.34 0.34 

Radio galaxies 1,542 0.04 0.48 

Total 10,000 0.58 0.55 



ergy, of sources classified as FSRQs and BL Lacs in 
our simulation with those of blazars included the radio 
sample of [2(|, which is the sample with the best de- 
termination of fpcak values currently available. The 
agreement is clearly quite good and reproduces well 
the fact that BL Lacs tend to have ^ pca k values sig- 
nificantly higher than FSRQs. 

4.2. X-ray flux limited survey 

Table ITT1 summarizes our main results. About 2/3 
of our blazars are classified as BL Lacs, which is con- 
sistent with the value of 731^% in the EMSS sample. 
As in the radio case, a small fraction (15%) of the sim- 
ulated blazars are misclassified as radio galaxies. The 
mean redshifts for our simulated FSRQs and BL Lacs 
are in reasonable agreement with the EMSS blazar 
sample values ~ 1 and ~ 0.37. We note that, unlike 
the WMAP5 sample, the EMSS sample is relatively 
small (56 sources) and therefore a detailed compari- 
son is hampered by the small number statistics. 79% 



of our BL Lacs have a redshift determination, in good 
agreement with the EMSS value of 93±g%. Although 
we assumed that all non-evolving sources do not have 
a standard accretion disk, 30% of the BL Lacs pos- 
sess one and are classified as BL Lacs only because 
their emission lines are swamped by the non-thermal 
continuum. The smaller fraction of X-ray selected BL 
Lacs with disks in our simulations, as compared to 
radio-selected ones, is in accordance with the fact that 
fewer EMSS BL Lacs have emission lines clearly de- 
tectable in their optical spectra than, for example, 
1 Jy BL Lacs [53, H . As regards (V/V m ), our 
simulated mean values agree reasonably well with the 
EMSS ones of 0.67 ± 0.08 and 0.42 ± 0.05 for FSRQs 
and BL Lacs respectively, derived using the samples 
described in [33[ . Fig. |4]compares the distributions of 
t'pcak of FSRQs and BL Lacs in our simulation with 
those of blazars belonging to the soft X-ray sample 
of [13], which includes Planck, Swift and Fermi ob- 
served blazars and it is therefore probably the sample 
with the best determination of fpeak values currently 
available. Our simulations reproduce the fact that BL 
Lacs have much higher z^ pca k values than FSRQs. 

4.3. A blazar sequence? 

The existence of an anti-correlation between bolo- 
metric luminosity and ^ pca k, the "blazar sequence" , 
has been discussed since first proposed by [1J| and 
[e.g., M M M M M O- Figure [5] shows 
our radio and X-ray selected simulated blazars in the 
log(^p Ca k) - log(^Ly(5 GHz)) plane. This reproduces 
the plot used by to propose the existence of the 
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Figure 5: The FSRQs (open squares) and BL Lacs (circles) of our radio flux density limited (/ r > 0.9 Jy, red symbols) 
and X-ray flux limited (/ x > 5 x 10~ 13 erg cm -2 s" 1 . green symbols) simulated samples, plotted in the log(^)-log(z/L„) 
plane, the portion of the parameter space used by [TI| to show the existence of the "blazar sequence" by comparing 
radio (/ r > 1 — 2 Jy) and X-ray selected (/ x Si 10 _12 erg cm -2 s _1 ) complete samples. Filled circles represent BL Lacs 
with very weak lines (EW < 2 A) or completely featureless which, in a real survey would most probably not have a 
measured redshift, and therefore would not appear in the plot. 



blazar sequence based on the correlation shown in 
this plane by FSRQs and BL Lacs discovered in shal- 
low radio surveys (2 and 1 Jy samples) and BL Lacs 
found in the X-ray flux limited Einstein slew survey 
{fx ~ 10 _12 erg cm~ 2 s _1 ). Indeed, considered to- 
gether the simulated radio and X-ray selected blazars 
display a broad correlation with radio selected FS- 
RQs and BL Lacs (red open squares and open circles) 
mostly filling the top left and central part of the dia- 
gram and X-ray selected BL Lacs (green open circles) 
mostly confined to the lower right corner of the plot. 
This particular positioning of the points (bright FS- 
RQs of the LSP type vs. fainter HSP BL Lacs) is 
not due to any intrinsic correlation between luminos- 
ity and fpeak but results from the fact that bright 
radio sources are mostly drawn from the high end of 
the blazar luminosity function, while BL Lacs in X- 
ray flux limited samples are mostly high v pea k sources 
(intrinsically rare) drawn from the low end of the lu- 
minosity function where the source density is largest. 
The most important difference between Fig. [5] and 
the diagram of fill ] is in the high-luminosity - high 
t'peak part, where most of the radio and X-ray selected 
objects with no redshift (red and green filled points) 
are located. These sources could not be plotted by 



since the luminosity of blazars without redshift 
cannot be estimated. This made the top right part of 
the (Tl| diagram empty, thus contributing to making 
the data points look like a sequence. 



5. DISCUSSION 
5.1. Evolution 

A long-standing blazar puzzle is the difference in 
redshift distribution and cosmological evolution be- 
tween FSRQs and BL Lacs, selected both in the radio 
and in the X-ray band. BL Lacs are mostly located at 
low redshifts and exhibit moderate, or even negative, 
evolution, while FSRQs evolve strongly just like ra- 
dio quiet QSOs and show a redshift distribution that 
peaks at z > 1 0, [H, HI Hi], H ■ Our simulations 
reproduce quite well both of these findings (see Fig. 
[2] and Tables Q] and |TT| implying that they are due 
to heavy selection effects. Most of the simulated BL 
Lacs found in radio surveys are luminous objects with 
broad lines that are diluted by non-thermal radiation 
beyond the 5 A EW limit (many of them just below, 
thus allowing a measurement of their redshift), while 
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the BL Lacs found in simulated X-ray surveys typi- 
cally show high i^p Ca k values (and therefore are X-ray 
bright) and are drawn from the low-power end of the 
radio luminosity function where non-evolving FR Is 
are preferentially found (Section 13. If) . We note that 
the (V/V m ) of radio selected BL Lacs is not too dif- 
ferent from that of FSRQs, while the (V/V m ) of X- 
ray selected BL Lacs is significantly lower, as found 
in real surveys ( Section l4~Tj) . Another interesting out- 
come of our simulations is the fact that X-ray selected 
BL Lacs with progressively larger values of f pea k are 
characterized by lower and lower values of (V/Vm) (see 
Tab. HI]) . This is in full agreement with the puzzling, 
and so far unexplained, results of [17( and [37| who 
reported that the (V/V m ) of BL Lacs is a function of 
their X-ray-to-radio flux ratio (which in turn depends 

On I/pcak)- 




1 2 3 

Redshift of featureless BL Lacs 

Figure 6: The redshift distribution of the BL Lacs that 
show a featureless spectrum in our simulation of a radio 
flux density limited survey and that in a real survey 
would have no redshift determination. 



5.2. is pcak distribution 

Recent results, based on radio and 7-ray surveys, 
have revealed that BL Lacs, on average, display a dis- 
tribution of ^p C ak energies, which is shifted to values 
higher than those of FSRQs 0, H(| , expanding on the 
well-known fact that high ^ poa k objects (HSPs) are 
always BL Lacs. This experimental difference is well 
reproduced in our simulations (see Fig. [3| which give 
(log(^poak)) - 12.9 for FSRQs and (log(iy pcak )) = 14.1 
for BL Lacs for the case of a radio survey. This dis- 
tinction is due to the fact that blazars with higher 
^peak values produce more non-thermal optical light 
than low f pea k sources, diluting more easily the broad 
line component, and are therefore classified more fre- 
quently as BL Lac objects. We note that this has 
been interpreted in the literature as an intrinsic phys- 
ical difference between LSP (detected mostly in the 
radio band) and HSP (detected mostly in the X-ray 
and 7-ray band) blazars due to the fact that HSPs are 
observationally characterized by a low intrinsic power 
and external radiation field, given their very weak or 
absent emission lines. As a consequence, cooling was 
thought to be less dramatic in HSP allowing particles 
to reach energies high enough to produce synchrotron 
emission well into the X-ray band [13j . In our scenario, 
instead, all sources have exactly the same chance of 
being HSP or LSP (that is, the value of 7 pea k is drawn 
from a distribution independently of luminosity) and 
the very different z^ poa k distributions observed in ra- 
dio and X-ray surveys arise from the strong selection 
effect discussed in Section [4731 and the emission line di- 
lution mentioned above. Our simulations predict the 
existence of a significant number of BL Lacs with red- 
shift that cannot be measured, which can occur when 
both iSp C ak and radio power are so large that dilution 
becomes extreme. For example, ~ 81% of our sim- 
ulated sources with ^ pca k > 10 15 Hz and P r > 10 26 
W/Hz have no redshift. This is consistent with the 



fact that most BL Lacs in current 7-ray selected sam- 
ples have no measured redshift, as Fermi is known to 
preferentially select high j/ pea k BL Lacs P, H[ . This ef- 
fect is also clearly shown in Fig. [5] where most of the 
simulated blazars with no measurable redshift (filled 
circles with light colors) occupy the top right part of 
the diagram. Fig. [5] shows the intrinsic redshift dis- 
tribution of these featureless BL Lacs for the case of 
our simulation of a radio flux density limited survey. 

5.3. What is a BL Lac? 

It was originally suggested that the absence of broad 
lines in BL Lacs was due to a very strong, Doppler- 
boosted synchrotron continuum 6]. In the years fol- 
lowing that paper observations of various BL Lacs, 
mostly selected in the X-ray band, showed that in 
many cases their optical spectrum was not swamped 
by a non-thermal component, as host galaxy features 
were very visible, and it was th oug ht that most BL 
Lacs had intrinsically weak lines [45| . We have shown 
here that these two possibilities are not mutually ex- 
clusive and indeed arc both viable, depending on radio 
power and, therefore, on the band of selection. One 
important consequence of our scenario is that objects 
so far classified as BL Lacs on the basis of their ob- 
served weak, or undetectable, emission lines belong 
to two physically different classes: intrinsically weak- 
lined objects, whose parents are LERGs/FR Is (more 
common in X-ray selected samples, since they reach 
lower radio powers) and heavily diluted broad-lined 
sources, which are beamed HERGs/FR lis (more fre- 
quent in radio selected samples). Therefore, while the 
non-thermal engine is probably the same, the ther- 
mal one is obviously different. This solves at once the 
issue of the FSRQ/BL Lac transition objects and of 
the many differences between BL Lacs selected in the 
radio and X-ray bands, which include line strength, 
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extended radio emission and morphology, and evolu- 
tion [e.g. l38l. and references therein]. It also implies 
that BL Lacertae, the prototype of the class, is ac- 
tually not a BL Lac but an FSRQ with its emission 
lines swamped by the jet. This hypothesis, which ex- 
plains also many other open issues of blazar research, 
is testable. Our simulations, in fact, imply that the 
majority of sources in high-flux density radio-selected 
samples are identified as BL Lacs only because all 
lines have EW < 5 A in the optical observing win- 
dow. These sources should show a strong (EW > 5 
A) Ha line in the near or mid-infrared, since this is 
the strongest emission line, and would then be con- 
sidered FSRQs. HERG/FR II BL Lacs, then, should 
be easily recognizable through infrared spectroscopy. 
These sources, as expected, are also more dominant at 
higher radio powers: for example, according to our ra- 
dio simulation the fraction of BL Lacs with standard 
accretion disks is only ~ 3% for P v < 10 26 W/Hz but 
becomes ~ 98% above this value. This also means 
that fainter radio-selected samples of BL Lacs should 
be more and more similar to X-ray selected ones, apart 
from their t'peak values, since higher ^ pca k are easier to 
detect in the X-rays [30] ■ The implications of our hy- 
pothesis for unified schemes is quite straightforward: 
the parent population of BL Lacs need to include both 
LERGs/FR Is and HERGs/FR lis, while that of FS- 
RQs is made up of HERGs/FR lis only. This should 
have only a small effect, for example, on the LF fit- 
ting done by [H, as HERG/FR II BL Lacs make up 
the high-power end of the radio LF while most of the 
number density comes from LERG/FR I BL Lacs. 

5.4. Blazar classification 

Our new scenario has strong implications on blazar 
classification. If the relevant physical distinction for 
radio sources is between LERGs (mostly FR Is) and 
HERGs (FR lis), for most purposes then (LFs, evo- 
lution, etc.) HERG/FR II BL Lacs should be simply 
grouped with FSRQs. How does one distinguish in 
practice HERG/FR II BL Lacs from LERG/FR I BL 
Lacs? This is simple in the presence of any (even 
weak) broad lines or for transition objects. In other 
cases (e.g. completely featureless spectrum or in pres- 
ence of absorption features) there is no easy way to 
distinguish between the two subclasses although, for 
example, radio power and/or morphology could help. 
Given the paucity of known LERG/FR Is at rela- 
tive high redshifts, X-ray selected (and also fainter 
radio-selected) BL Lac samples are also useful in se- 
lecting such sources, which are very relevant also for 
the study of the so-called "AGN feedback" and the 
role that AGN radio emission plays in galaxy evo- 
lution through the so-called "radio-mode" accretion 
Q. It should also be clear that BL Lacs can be 
used to study the broader issue of the relationship 



between LERGs and HERGs, including their evolu- 
tion. A by-product of our simulations has also been 
the realization that some sources classified as radio- 
galaxies do not have their jets oriented at large angles 
with respect to the line of sight, as expected, but are 
instead moderately beamed blazars with their non- 
thermal emission swamped by the galaxy (note that 
none of these objects has a standard accretion disk). 
These sources, which are all local (~ 90% at z < 0.07) 
should be recognizable by their blazar-like SEDs and 
indeed some of them have already been identified by 
H, [3 and 0- Recently, a new classification scheme 
has been proposed to divide BL Lacs from FSRQs, 
which is based on the broad line region (BLR) lumi- 
nosity in Eddington units and set at a dividing value of 
iBLR/iEdd ~5x 10~ 4 15]. This turns out to be also 
the value, which separates radiativcly efficient (i.e., 
standard accretion disks) from radiatively inefficient 
(i.e., ADAFs) regimes, and therefore coincides with 
our HERG/FR II - LERG/FR I division. Therefore, 
[l5l ] are also suggesting that HERG/FR II BL Lacs 
belong with the FSRQs. 



6. TESTS AND PROSPECTS 

Very recently, redshift constraints for 103 blazars 
from the Fermi 2LAC catalogue [3j have been 
derived by fitting SED templates to their UV- 
to-near-IR multi-band photometry obtained quasi- 
simultaneously with Swift/UVOT and GROND jH]. 
This was done using the attenuation due to neutral 
hydrogen along the line of sight at the Lyman limit to 
estimate the redshift of the absorber. Eleven of these 
objects have z p hot > 1-2. Some of us have studied the 
SEDs of these sources using quasi-simultaneous near- 
IR to X-ray data. Four blazars turned out to be of 
a type never seen before but with properties we were 
expecting for some of the BL Lacs without redshift: 
large ^ poa k (~5x 10 15 Hz) and high-power (see filled 
circles in Fig. [S]) [12]. Given their featureless optical 
spectra, these sources are therefore most likely high- 
redshift FSRQs with their emission lines swamped by 
the jet, as predicted by our hypothesis. 

In this paper we limited our simulations to the ra- 
dio and the X-ray bands where SSC is a fair approx- 
imation of the observed non-thermal emission. The 
properties of 7-ray detected blazars are instead not 
consistent with simple SSC models [e.g.Q, and almost 
half of the radio and X-ray selected LSP blazars are 7- 
ray quiet [e.g. 20]. The present approach must there- 
fore be integrated with additional information about 
the properties of the inverse Compton emission be- 
fore it can be used to simulate 7-ray surveys. We are 
planning to extend our simulations to the 7-ray band 
by taking into account the recent results of [2(j] who 
determined the 7-ray properties of blazar samples se- 
lected in different bands. 
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